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Temperature  and  humidity  fluctuations  in  operating  fuel  cells  impose  significant  biaxial  stresses  in  the 
constrained  proton  exchange  membranes  (PEMs)  of  a  fuel  cell  stack.  The  strength  of  the  PEM,  and  its 
ability  to  withstand  cyclic  environment-induced  stresses,  plays  an  important  role  in  membrane  integrity 
and  consequently,  fuel  cell  durability.  In  this  study,  a  pressure  loaded  blister  test  is  used  to  characterize 
the  biaxial  strength  of  Gore-Select®  series  57  over  a  range  of  times  and  temperatures.  Hencky’s  classi¬ 
cal  solution  for  a  pressurized  circular  membrane  is  used  to  estimate  biaxial  strength  values  from  burst 
pressure  measurements.  A  hereditary  integral  is  employed  to  construct  the  linear  viscoelastic  analog  to 
Hencky’s  linear  elastic  exact  solution.  Biaxial  strength  master  curves  are  constructed  using  traditional 
time-temperature  superposition  principle  techniques  and  the  associated  temperature  shift  factors  show 
good  agreement  with  shift  factors  obtained  from  constitutive  (stress  relaxation)  and  fracture  (knife  slit) 
tests  of  the  material. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Strength  characterization  of  proton  exchange  membranes 
(PEMs)  is  critical  for  lifetime  prediction  of  the  membranes  in 
fuel  cell  applications.  The  Department  of  Energy  has  set  durabil¬ 
ity  goals  for  PEM  transportation  applications  at  a  5000  h  service 
life  by  2015  [1].  These  goals  would  hold  PEM  fuel  cell  vehicles  to 
the  same  service  life  standards  of  today’s  current  internal  com¬ 
bustion  engine  systems.  However,  current  fuel  cell  demonstration 
vehicles  have  yet  to  reach  the  DOE  target  [2].  One  reason  is  that 
under  operating  conditions  and  due  to  non-continuous  usage,  the 
fuel  cell  environment  experiences  both  temperature  and  humid¬ 
ity  fluctuations.  When  exposed  to  such  hygrothermal  changes,  the 
membrane,  which  is  clamped  along  the  lands  and  seal  gasket  of 
fuel  cells  and  constrained  by  the  catalyst  and  gas  diffusion  layers,  is 
expected  to  be  under  biaxial  stresses.  Lai  et  al.  has  shown  that  these 
imposed  cyclic  stresses  can  initiate  crack  formation,  which  leads 
to  gas  crossover  and  subsequent  fuel  cell  failure  [3].  In  a  previous 
work,  Dillard  et  al.  [4]  showed  that  pressure  loaded  blister  testing 
is  a  viable  method  to  impose  these  biaxial  stress  states  experimen¬ 
tally  and  Li  et  al.  reported  the  use  of  pressure-loaded  blisters  to 
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characterize  the  fatigue  and  creep  to  leak  behaviors  of  PEMs  [5].  In 
particular,  pressure-ramped  blister  testing  provides  biaxial  leak¬ 
age  or  burst  strength  estimates  which  could  be  more  descriptive  of 
failure  modes  seen  in  working  fuel  cells  than  traditional  strengths 
obtained  by  uniaxial  tensile  testing.  Since  biaxial  burst  strength 
values  could  provide  insight  into  evaluating  membrane  durability, 
a  means  of  experimentally  and  analytically  obtaining  these  values 
while  incorporating  accurate  material  behavior  is  desired. 

In  early  work  with  pressurized  circular  membrane  testing, 
Hencky  showed  that  stresses  within  the  pressurized  blister  could 
be  exactly  determined  for  a  linear  elastic,  homogeneous,  isotropic 
membrane,  using  a  power  series  solution  [6].  The  radial  and  tan¬ 
gential  stresses,  or  and  respectively,  are  given  [6,7]  as: 
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where  E  is  the  elastic  modulus  of  the  material,  p  is  the  applied  pres¬ 
sure,  a  is  the  free  radius  of  the  blister,  h  is  the  thickness  of  the 
blister,  and  B2/<  is  the  coefficients  to  be  determined  through  applica¬ 
tion  of  boundary  conditions.  By  assuming  zero  residual  mechanical 
strain,  a  reasonable  approximation  of  or  and  ae  at  the  center  of  a 
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pressurized  PEM  blister  test  is  given  by, 


Or  =Oq=  — 
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where  1.777  for  an  assumed  Poisson’s  ratio  value  of  v=0.4. 
Since  Hencky’s  solution  assumes  a  linear  elastic  material,  modifi¬ 
cations  are  needed  to  address  the  viscoelastic  nature  of  polymers 
like  PEM.  To  accurately  explore  time  and  temperature  dependence 
of  biaxial  strength,  material  behavior  due  to  environmental  condi¬ 
tions  must  also  be  considered. 


2.  Theory 

The  viscoelastic  response  of  polymers  such  as  PEMs  exhibits  a 
dependence  on  temperature,  time,  and  humidity.  Early  work  with 
polymers  led  to  the  introduction  of  the  time-temperature  superpo¬ 
sition  principle  (TTSP),  which  allows  for  the  transformation  of  data 
from  one  temperature  to  another  via  horizontal  shifts  on  a  loga¬ 
rithmic  time  axis  [8].  Alternately,  short  term  viscoelastic  properties 
such  as  relaxation  modulus  collected  at  different  temperatures  can 
be  superimposed  into  a  master  curve  describing  behavior  at  times 
both  shorter  and  longer  than  experimentally  accessible.  However, 
as  noted  in  [8],  the  simple  TTSP  is  not  applicable  to  all  polymers. 
Relevant  to  this  study,  success  has  been  demonstrated  in  applying 
TTSP  to  Gore-Select®  series  57.  Specifically,  Patankar  et  al.  [9]  stud¬ 
ied  the  effects  of  both  temperature  and  humidity  on  Gore-Select® 
series  57  and  presented  an  initial  linear  viscoelastic  characteriza¬ 
tion  in  the  form  of  a  relaxation  modulus  master  curve  along  with 
hygrothermal  shift  factors. 

The  behavior  of  these  viscoelastic  membranes  as  relevant  to  our 
study  is  illustrated  in  Fig.  1,  as  a  plot  of  the  Prony  series  approx¬ 
imation  of  the  relaxation  modulus  master  curve  obtained  from 
stress  relaxation  tests  conducted  on  the  Gore-Select®  series  57 
membranes  tested  in  a  dynamic  mechanical  analyzer  (DMA)  as 
presented  in  [9].  It  is  important  to  note  that  the  data  of  Fig.  1 
was  collected  under  dry  conditions  to  best  simulate  the  uncon¬ 
trolled  humidity  of  the  high  temperature  blister  tests  conducted 
above  and  is  not  identical  to  the  curve  presented  in  [9].  The  regions 
denoted  on  the  plot  show  the  modulus  variations  within  the  time 
and  temperature  ranges  of  interest  to  this  study. 

It  has  been  shown  that  the  elastic-viscoelastic  correspon¬ 
dence  principle  can  be  used  to  derive  linear  viscoelastic  solutions 
from  known  linear  elastic  solutions  provided  the  geometries  and 


boundary  conditions  are  identical  [10].  Also  known  as  Alfrey’s 
correspondence  principle  the  process  of  obtaining  these  linear  vis¬ 
coelastic  solutions  is  described  in  detail  in  [11].  Each  traction  or 
displacement  variable  appearing  in  the  analogous  elastic  solution 
is  replaced  by  the  Laplace  transform  of  that  variable.  Next,  any  elas¬ 
tic  constant  appearing  in  the  solution  is  replaced  by  the  product 
of  the  Laplace  variable  s  and  the  Laplace  transform  of  the  time 
dependent  analog  to  the  elastic  constant.  Having  now  obtained 
the  linear  viscoelastic  solution  in  the  Laplace  domain,  we  then 
take  the  inverse  transform  to  complete  the  procedure  and  pro¬ 
vide  the  analogous,  time  dependent,  linear  viscoelastic  solution.  In 
application  to  Hencky’s  solution  described  above,  the  only  variable 
requiring  transformation  into  the  Laplace  domain  is  the  applied 
pressure.  Similarly,  the  only  elastic  constant  appearing  directly  in 
Hencky’s  solution  is  the  Young’s  modulus,  E,  whose  corresponding 
time  dependent  viscoelastic  property  is  the  stress  relaxation  modu¬ 
lus  E(t).  Although  this  principle  is  typically  presented  in  traditional 
textbook  examples,  such  as  simple  shaft  torsion,  this  methodology 
has  been  shown  broadly  useful  in  solving  more  complex  mechanics 
problems,  such  as  in  the  indentation  field  [12,13]. 

In  applying  the  correspondence  principle  to  Eq.  (2)  and  taking 
the  inverse  Laplace  transform,  one  obtains  a  hereditary  integral  for¬ 
mulation  in  which  the  stress  relaxation  modulus,  E(t),  is  convoluted 
with  the  pressure  input  as  the  only  other  time  dependent  parame¬ 
ter.  These  modifications  give  ar  and  ae  at  the  center  of  a  pressurized 
PEM  blister  as 


ar(t)  =  ae(t)=§.(£)2/3^  (£(f-j9))1/3.d[p2d/3(/i)1^  (3) 

Building  upon  the  above  concepts,  the  objective  of  this  study  is  to 
characterize  the  biaxial  burst  strength  of  Gore-Select®  series  57 
PEM  over  a  range  of  time  and  temperature  values  using  a  pres¬ 
surized  blister  test.  A  hereditary  integral  formulation  of  Hencky’s 
classical  solution  will  be  employed  using  stress  relaxation  data 
obtained  from  the  methodology  in  [9].  The  TTSP  will  be  inde¬ 
pendently  applied  to  the  collected  data  and  the  temperature  shift 
factors  obtained  will  be  compared  to  the  shift  factors  collected  from 
constitutive  (stress  relaxation)  and  fracture  (knife  slitting)  tests  as 
described  in  [9,14].  From  application  of  TTSP,  master  curves  will 
be  generated  that  describe  the  effects  of  time  and  temperature  on 
biaxial  burst  strength. 

3.  Experimental  procedure  and  analysis 


Time,  sec 

Fig.  1.  Plot  of  Prony  series  approximation  of  the  relaxation  modulus  obtained  from 
DMA  tests  of  Gore-Select®  series  57  under  dry  conditions  with  temperature  regions 
denoted  as  applicable  to  the  time  range  of  the  data  utilized  in  this  study. 


The  pressurized  blister  testing  method  employed  for  PEMs  in 
previous  work  [4]  was  expanded  in  this  study  to  incorporate  a 
range  of  test  temperatures  in  order  to  investigate  the  applicabil¬ 
ity  of  the  TTSP  in  pressure  loaded  blister  tests  over  temperatures 
relevant  to  operating  fuel  cells.  Specifically,  tests  were  conducted 
on  Gore-Select®  series  57  (Gore  Fuel  Cell  Technology,  Elkton, 
MD)  at  temperatures  of  70,  80,  and  90  °C  and  readings  of  the 
pressure  at  burst  and  elapsed  time  at  burst  were  collected.  The 
experimental  procedure,  including  sample  preparation,  heating 
to  the  desired  test  temperature,  and  pressurizing  until  bursting 
occurs,  was  described  in  detail  in  [4]  and  is  only  briefly  reiter¬ 
ated  here.  Gore-Select®  series  57  is  an  18fxm-thick  version  of 
the  perfluorosulfonic  acid  (PFSA)  proton  exchange  membrane  with 
a  layer  of  micro-reinforced,  composite  expanded  polytetrafluo- 
roethylene  (ePTFE).  These  micro-reinforced  membranes  have  been 
shown  to  have  improved  tear  resistance  and  dimensional  stability 
compared  to  homogeneous  cast  membranes  under  some  condi¬ 
tions  [5,15,16].  Gore-Select®  series  57  membrane,  as  received,  is 
sandwiched  between  two  backings.  One  backing  was  removed 
before  specimen  preparation,  while  the  other  backing  was  retained 
throughout  specimen  fabrication  to  form  consistent  samples.  Circu- 
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lar  specimens  were  constructed  using  the  front  ferrules  of  standard 
Swagelok®  tube  fittings.  A  thin  coat  of  Devcon®  two-part  5-min 
epoxy  was  applied  on  the  rim  of  the  front  ferrules,  which  were 
then  inverted  and  placed  onto  a  sheet  of  the  membrane.  After  cur¬ 
ing  at  room  temperature,  individual  specimens  were  cut  out  using 
a  razor  blade,  the  second  backing  was  removed,  and  the  specimens 
were  placed  in  a  desiccator  until  specimen  testing.  To  mount  a 
specimen  in  preparation  for  testing,  the  front  ferrule  supporting 
a  just  taut  membrane  was  placed  into  a  Swagelok®  ferrule  nut  and 
the  assembly  was  screwed  onto  a  Swagelok®  reducing  union.  Once 
mounted  on  the  pressure  line,  specimens  were  placed  in  an  oven  for 
environmental  control  and  were  pressurized  using  a  KDScientifc® 
Series  230  syringe  pump  with  atmospheric  air  as  the  working  fluid. 
A  schematic  of  the  experimental  setup  is  shown  in  Fig.  2a  and  a 
photograph  of  the  test  in  progress  is  shown  in  Fig.  2b,  respectively. 

Raw  pressure  versus  time  data  were  collected  from  each  test, 
where  test  initiation  was  marked  by  the  first  noticeable  increase  in 
pressure  reading  and  bursting  was  denoted  by  catastrophic  pres¬ 
sure  loss.  Successful  testing  of  catastrophic  burst  failures  was  also 
verified  in  a  qualitative  way  by  the  audible  popping  noise  that 
accompanied  such  failures.  Raw  data  from  a  representative  sample 
is  plotted  as  pressure  versus  time  in  Fig.  3. 

From  the  collected  readings  from  similar  tests,  burst  pressures 
and  times  to  failure  were  obtained  for  a  series  of  syringe  displace¬ 
ment  rates  and  test  temperatures.  For  the  analysis  of  this  data,  Eq. 
(3)  was  utilized  extensively  using  temperature  appropriate  relax¬ 
ation  modulus  data  given  by  the  Prony  series  approximation  for 
Gore-Select®  series  57  presented  in  Fig.  1.  Although  this  approach 
proves  to  be  a  vast  improvement  over  the  assumptions  of  linear 
elastic  behavior  within  the  classical  Flencky  solution,  it  is  imper¬ 
ative  to  objectively  recall  assumptions  made  in  formulating  this 


Fig.  2.  Illustrations  of  (a)  experimental  schematic  diagram  and  (b)  blister  test  in 
progress  [4]. 


Time,  sec 


Fig.  3.  Raw  pressure  versus  time  data  from  representative  sample  of  a  Gore-Select® 
series  57  blister  test  conducted  at  80  °C  with  a  syringe  pump  infusion  rate  of 
12  mL  min-1.  Associated  time  and  pressure  at  burst  are  illustrated. 


analogous  linear  viscoelastic  solution.  In  the  procedure,  the  relax¬ 
ation  modulus  E(t )  was  the  only  time  dependent  material  property 
used  in  the  modification  and  was  convoluted  in  the  integral  with  the 
applied  pressure  as  a  function  of  time,  p(t).  Admittedly,  the  thick¬ 
ness  will  also  be  a  time  dependent  parameter  as  the  membrane 
continues  thinning  throughout  the  duration  of  the  test.  However, 
appropriate  treatment  of  the  thickness  becomes  quite  complicated 
as  it  would  involve  not  only  applied  stress  levels  at  each  instant  in 
time  but  also  a  viscoelastic  characterization  of  Poisson’s  ratio,  v(t ), 
which  is  not  well  known  or  published  for  Gore-Select®  series  57. 
To  significantly  reduce  computational  efforts,  a  constant  value  for 
the  thickness  of  the  membrane  was  assumed. 

The  motivation  behind  developing  and  employing  this  hered¬ 
itary  integral  formulation  was  not  just  based  on  a  desire  to 
incorporate  the  time  dependent  relaxation  modulus.  The  primary 
reason  resulted  from  anomalous  results  found  from  an  attempted 
application  of  TTSP  to  data  obtained  with  an  assumed  constant  elas¬ 
tic  modulus.  The  results  were  temperature  shift  factors  that  were 
surprisingly  large,  approximately  twice  the  values  obtained  from 
DMA  testing.  It  was  known  that  as  the  hereditary  integral  formula¬ 
tion  accounts  for  time  dependency  of  the  relaxation  modulus,  the 
overall  slope  of  the  biaxial  burst  strength  versus  time  curve  would 
become  larger  in  magnitude  compared  to  data  with  assumed  con¬ 
stant  elastic  modulus.  In  application  of  TTSP,  the  steeper  slopes  of 
the  hereditary  integral  method  required  less  horizontal  shifting  to 
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Fig.  4.  Comparison  of  slopes  in  curves  obtained  via  constant  elastic  modulus  and 
Hencky  with  hereditary  integral  methods. 
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Fig.  5.  Burst  strength  data,  calculated  via  Hencky’s  solution  with  a  hereditary 
integral,  shifted  to  form  a  master  curve  of  Gore-Select®  series  57  obtained  using 
traditional  shifting  methods  and  referenced  at  70  °C. 

obtain  overlap,  thereby  reducing  the  size  of  the  temperature  shift 
factors.  This  key  concept  of  shift  factor  difference  due  to  general 
slope  difference  is  presented  visually  in  Fig.  4. 

The  data  reduction  process  outlined  above  provided  a  set  of 
burst  strength  versus  time  to  failure  curves  at  selected  tempera¬ 
tures.  Traditional  shifting  methods  were  utilized  to  form  a  burst 
strength  master  curve  with  a  reference  temperature  of  70  °C.  Tem¬ 
perature  shift  factors  generated  were  then  compared  to  shift  factors 
obtained  by  stress  relaxation  and  published  knife  slit  tests  in  [14] 
under  dry  conditions  to  evaluate  consistency. 

4.  Results  and  discussion 


Fig.  7.  Plot  of  temperature  shift  factors  obtained  using  methodology  from  [9,14]  for 
stress  relaxation  and  knife  slit  tests  and  compared  with  shift  factors  obtained  in  this 
study. 


Table  1 

Table  of  temperature  shift  factors  obtained  using  methodology  from  [8,13]  for  stress 
relaxation  and  knife  slit  tests  and  compared  with  shift  factors  obtained  in  this  study. 


70  °C 

80  °C 

90  °C 

Stress  relaxation  tests 

(Reference) 

-1.00 

-2.50 

Knife  slit  tests 

(Reference) 

-1.20 

-2.40 

Burst  strength 
tests— Hencky 

W/hereditary  integral 

(Reference) 

-1.00 

-2.05 

Burst  strength 
tests— Hencky 

W/constant  modulus 

(Reference) 

-1.50 

-3.10 

Application  of  the  hereditary  integral  to  Hencky’s  solution  led  to 
a  collection  of  burst  strength  values  and  associated  times  to  failure 
over  the  temperature  regions  of  interest.  This  data  was  then  shifted 
using  the  traditional  method  of  choosing  a  reference  data  set  and 
sliding  subsequent  data  sets  horizontally  on  a  logarithmic  time  axis 
until  visual  alignment  with  the  previous  data  set  was  obtained.  No 
vertical  shifts  were  employed.  This  method  of  visual  shifting  was 
used  with  a  chosen  reference  temperature  of  70  °C  to  form  the  burst 
strength  master  curve  featured  in  Fig.  5. 

As  mentioned  previously,  the  use  of  a  hereditary  integral  with 
Hencky’s  solution  gives  an  exact  linear  viscoelastic  solution  and  is 


Time,  sec 

Fig.  6.  Comparison  plot  of  burst  strength  master  curves  based  on  both  Hencky’s 
solution  with  a  hereditary  integral  and  constant  modulus  stress  analysis  methods. 


considered  an  improvement  over  the  constant  modulus  approach. 
For  completeness,  a  comparison  plot  of  the  burst  strength  mas¬ 
ter  curves  obtained  via  hereditary  integral  formulation  of  Hencky’s 
solution  and  constant  modulus  approach  is  shown  in  Fig.  6. 

One  criterion  for  evaluating  the  appropriateness  of  TTSP  is  the 
consistency  of  shift  factors  obtained  for  different  material  proper¬ 
ties  [17].  It  is  therefore  encouraging  to  note  that  in  comparisons 
of  these  results  with  published  shift  factors  obtained  from  stress 
relaxation  tests  and  knife  slitting  tests,  the  Hencky  solution  with 
hereditary  integral  was  found  to  be  in  reasonable  agreement  [9,14]. 
This  suggests  that  the  same  molecular  relaxation  processes  that 
govern  stress  relaxation  behavior  and  fracture  in  membranes  also 
control  the  time  dependent  failures  during  biaxial  loading.  Fig.  7 
and  Table  1  illustrate,  in  graphical  and  tabular  forms,  the  agreement 
of  the  different  temperature  shift  factors  obtained  for  Gore-Select® 
series  57.  Development  of  a  burst  strength  master  curve  that  pro¬ 
vides  shift  factors  comparable  to  those  obtained  from  shifting 
constitutive  data  confirms  the  appropriateness  of  a  viscoelastic 
framework  for  characterizing  PEM  behavior  and  suggests  the  possi¬ 
bility  of  basing  lifetime  prediction  efforts  on  such  accelerated  tests. 
Although  shift  factor  agreement  does  not  guarantee  success,  it  does 
support  attempts  to  extend  strength-life  curve  results  from  short 
time  scales  into  longer  lifetime  investigations.  This  is  the  focus  of 
an  ongoing  study. 

5.  Conclusions 

Pressure  loaded  blister  tests  were  used  to  characterize  the 
biaxial  strength  of  commercially  available  Gore-Select®  series  57 
membrane  over  a  range  of  times  to  burst  and  temperatures  of 
interest.  The  time  temperature  superposition  principle  was  used  to 
obtain  biaxial  strength  master  curves  using  the  traditional  visual 
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shifting  method  on  data  reduced  via  the  Hencky  solution  with 
both  a  constant  modulus  and  a  hereditary  integral  linear  vis¬ 
coelastic  solution;  the  resulting  shift  factors  from  the  hereditary 
integral  method  were  quite  consistent  with  temperature  shift  fac¬ 
tors  obtained  from  stress  relaxation  and  knife  slitting  tests. 

While  considering  the  encouraging  insight  and  progress  gained 
from  this  study,  it  is  important  to  recall  that  the  data  reduction 
relies  heavily  on  the  linear  viscoelastic  material  characterization 
from  stress  relaxation  tests.  This  characterization,  which  provides 
no  insight  into  nonlinear  or  elastic-plastic  material  properties,  is 
still  utilized  in  this  study  to  treat  a  problem  that  likely  involves 
both  nonlinearity  and  plasticity.  Although  success  was  achieved 
via  temperature  shift  factor  agreement,  the  burst  strength  values 
presented  are  likely  lower  than  the  true  values  would  be  due  to 
localized  plasticity  and  the  material  nonlinearity.  In  spite  of  these 
approximations  and  omissions,  the  methodology  is  still  useful  as  a 
cost-  and  time-effective  tool  in  better  understanding  and  improv¬ 
ing  PEM  in  operating  fuel  cells.  In  particular,  these  tests  can  be 
used  to  quickly  screen  membranes  and  compare  the  biaxial  burst 
strength  of  one  material  to  another.  This  data  can  also  be  used 
to  investigate  the  critical  issue  of  membrane  durability  and  reli¬ 
ability.  Repeated  trials  at  identical  loading  rates  can  be  used  to 
begin  statistical  strength  studies  leading  to  material  reliability  con¬ 
stants  based  on  an  assumed  statistical  distribution.  Additionally, 
a  burst  strength  master  curve  under  ramped  pressure  conditions 
can  lend  insight  into  possible  relationships  between  applied  stress 
levels  and  subsequent  times  to  failure.  Given  a  suggested  relation¬ 
ship  for  time  to  failure  as  a  function  of  applied  stress,  a  damage 
evolution  law  can  be  applied  and  experimental  failures  under 
given  loading  conditions  can  be  predicted.  The  ability  to  develop 
these  techniques  in  such  a  simple  and  time-effective  test  would  be 
invaluable  if  proven  to  predict  more  complex  loading  histories.  If 
these  efforts  in  lifetime  prediction  model  are  pursued  and  demon¬ 
strated  with  success,  membrane  durability  could  be  investigated 
quickly  and  play  a  crucial  role  in  achieving  the  DOE  service  life 
targets  [18].  As  another  component  to  be  added  to  this  study,  a 
similar  characterization  could  be  performed  over  a  range  of  tem¬ 
peratures  and  relative  humidities  simultaneously  to  evaluate  the 
hygrothermal  dependence  of  the  biaxial  strength  of  Gore-Select® 
series  57.  Finally,  extending  all  of  the  applications  mentioned  into 
comparative  studies  of  other  commercially  available  membranes 
would  provide  valuable  insight  into  behavioral  differences  among 
the  products. 
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